
*** REPLICATION DATA FOR DILL, SAGAN, VALENTINO

*************************** FIGURE 2: nuke preference and nuke ethical by country and all conditions with conventional ethics (see below for Christian conditions)
graph drop _all
logit pref_nuke_A_bin treatment_A##country_ascending [pweight= weight] if experiment==1 & treatment_A!=5
set scheme s2mono
margins, over (treatment_A country_ascending)
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.8)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% prefer nuclear") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title("A. If you had to choose between one of the two military options described in the news article," "would you prefer the nuclear strike or the conventional strike?", size(small)) legend(off) xsize(17) ysize(11) fysize(51) name(scat1)
logit nuke_ethical_A_bin treatment_A##country_ascending [pweight= weight] if experiment==1 & treatment_A!=5
set scheme s2mono
margins, over (treatment_A country_ascending)
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.8)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% nuclear ethical") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title("B. Regardless of which option you preferred, how ethical or unethical" "do you think it would be if [country] decided to conduct the nuclear strike in this situation?", size(small)) legend(off) xsize(17) ysize(11) fysize(51) name(scat2)
logit conv_ethical_A_bin treatment_A##country_ascending [pweight= weight] if experiment==1 & treatment_A!=5
set scheme s2mono
margins, over (treatment_A country_ascending)
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.8)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% conv. ethical") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title("C. Regardless of which option you preferred, how ethical or unethical" "do you think it would be if [country] decided to conduct the conventional strike in this situation?", size(small)) legend(rows(1)) xsize(17) ysize(11) fysize(65) name(scat3)
graph combine scat1 scat2 scat3, imargin(vsmall) col(1)


** difference of means for both ethical and preference for figure 2
mean nuke_ethical_A_bin pref_nuke_A_bin, over(country treatment_A), [pweight=weight], coefleg
test _b[c.nuke_ethical_A_bin@2.country#1.treatment_A] = _b[c.nuke_ethical_A_bin@2.country#2.treatment_A]
test _b[c.nuke_ethical_A_bin@1.country#1.treatment_A] = _b[c.nuke_ethical_A_bin@1.country#2.treatment_A]
test _b[c.nuke_ethical_A_bin@3.country#1.treatment_A] = _b[c.nuke_ethical_A_bin@3.country#2.treatment_A]
test _b[c.nuke_ethical_A_bin@4.country#1.treatment_A] = _b[c.nuke_ethical_A_bin@4.country#2.treatment_A]

** pooling across conditions 1-4
mean pref_nuke_A_bin, over(country), [pweight=weight], coefleg, if treatment_A!=5


* interactions

* prefer nukes interactions
logit pref_nuke_A_bin treatment_A##country [pweight= weight] if treatment_A==1 | treatment_A==2
margins treatment_A#country
marginsplot, plotregion(margin(8 8 2 2)) plotopts(mlabel(_margin) mlabangle(45) yvarformat(%5.4f) mlabsize(vsmall)) xtitle("treatment group") xlabel(, angle(45) labsize(vsmall))
contrast treatment_A#rb2.country, effects
contrast treatment_A#rb3.country, effects
contrast treatment_A#rb4.country, effects

logit pref_nuke_A_bin treatment_A##country [pweight= weight] if treatment_A==2 | treatment_A==3
margins treatment_A#country
marginsplot, plotregion(margin(8 8 2 2)) plotopts(mlabel(_margin) mlabangle(45) yvarformat(%5.4f) mlabsize(vsmall)) xtitle("treatment group") xlabel(, angle(45) labsize(vsmall))
contrast treatment_A#rb2.country, effects
contrast treatment_A#rb3.country, effects
contrast treatment_A#rb4.country, effects

* nukes ethical country interactions (none found)
logit nuke_ethical_A_bin treatment_A##country [pweight= weight] if country!=5 &(treatment_A==1 | treatment_A==2) 
margins treatment_A#country
marginsplot, plotregion(margin(8 8 2 2)) plotopts(mlabel(_margin) mlabangle(45) yvarformat(%5.4f) mlabsize(vsmall)) xtitle("treatment group") xlabel(, angle(45) labsize(vsmall))
contrast treatment_A#rb2.country, effects
contrast treatment_A#rb3.country, effects
contrast treatment_A#rb4.country, effects

logit nuke_ethical_A_bin treatment_A##country [pweight= weight] if country!=5 &(treatment_A==2 | treatment_A==3) 
margins treatment_A#country
marginsplot, plotregion(margin(8 8 2 2)) plotopts(mlabel(_margin) mlabangle(45) yvarformat(%5.4f) mlabsize(vsmall)) xtitle("treatment group") xlabel(, angle(45) labsize(vsmall))
contrast treatment_A#rb2.country, effects
contrast treatment_A#rb3.country, effects
contrast treatment_A#rb4.country, effects

logit nuke_ethical_A_bin treatment_A##country [pweight= weight] if country!=5 &(treatment_A==1 | treatment_A==4) 
margins treatment_A#country
marginsplot, plotregion(margin(8 8 2 2)) plotopts(mlabel(_margin) mlabangle(45) yvarformat(%5.4f) mlabsize(vsmall)) xtitle("treatment group") xlabel(, angle(45) labsize(vsmall))
contrast treatment_A#rb2.country, effects
contrast treatment_A#rb3.country, effects
contrast treatment_A#rb4.country, effects

**conventional ethical interactions
logit conv_ethical_A_bin treatment_A##country_ascending male age_over60 relig_imp_bin college_grad support_conservative [pweight= weight] if (treatment_A==1 | treatment_A==2) & country_ascending!=5
margins treatment_A#country_ascending
marginsplot, plotregion(margin(8 8 2 2)) plotopts(mlabel(_margin) mlabangle(45) ylabel(#10,grid labsize(small)) yscale(range(0 1)) yvarformat(%5.2f) mlabsize(vsmall)) ytitle("% conventional ethical") xtitle("treatment group", size(3)) xlabel(, labsize(vsmall)) title("conventional ethical?", size(3)) legend(rows(1) size(2.8)) xsize(10) ysize(8)


*************************** FIGURE 3: Ranking of “Do Not Strike” Option, by Country and Condition CONDITIONS 1 and 4 ONLY

logit rank1_nostrike_A_bin treatment_A##country_ascending [pweight= weight] if experiment==1 & (treatment_A==1 | treatment_A==4)
set scheme s2mono
margins, over (treatment_A country_ascending)
mplotoffset, offset(0.35) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.8)) ylabel(#10,grid) yscale(range(0 1)titlegap(*12)) ytitle("% ranked no strike 1st") xtitle(" ") xlabel(, angle(0) labsize(small)) title("Rank the following options from 1 to 3," "where 1 indicates the option that you prefer the most" "and 3 indicates the option that you prefer the least.", size(small)) legend(rows(1)) xsize(11) ysize(9)

* compare means of 4-country average between conditon 1 and condition 4
mean rank1_nostrike_A_bin, over(country treatment_A), [pweight=weight], coefleg
test _b[c.rank1_nostrike_A_bin@5.country#1bn.treatment_A] = _b[c.rank1_nostrike_A_bin@5.country#4.treatment_A]


*************************** FIGURE 4: Preference for Nuclear Strike Option and Ethical Assessments by Compatriot Partiality and Death Penalty Support

graph drop _all
logit pref_nuke_A_bin treatment_A##own_vs_Libya_feeling_50_label male age_over60 relig_imp_bin college_grad support_conservative [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 & treatment_A!=5
set scheme s2mono
margins treatment_A#own_vs_Libya_feeling_50_label 
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.2)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% prefer nuclear") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title("If you had to choose between one of the two military options described in the news article," "would you prefer the nuclear strike or the conventional strike?" "     " "    "  "A. COMPATRIOT PARTIALITY", size(small)) legend(rows(1) size(small)) xsize(12) ysize(11) fysize(60) name(scat1)
logit pref_nuke_A_bin treatment_A##approve_deathpen male age_over60 relig_imp_bin college_grad support_conservative [pweight= weight] if experiment==1 & country!=5 & treatment_A!=5
set scheme s2mono
margins treatment_A#approve_deathpen
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.2)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% prefer nuclear") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title(" " "B. DEATH PENALTY", size(small)) legend(rows(1) size(small)) xsize(12) ysize(11) fysize(51) name(scat2)
graph combine scat1 scat2, imargin(medium) col(1)

**using "Be Like Country" measure of compatriot partiality
graph drop _all
logit pref_nuke_A_bin treatment_A##be_like_country_bin male age_over60 relig_imp_bin college_grad support_conservative [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 & treatment_A!=5
set scheme s2mono
margins treatment_A#be_like_country_bin 
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.2)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% prefer nuclear") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title("If you had to choose between one of the two military options described in the news article," "would you prefer the nuclear strike or the conventional strike?" "     " "    "  "A. COMPATRIOT PARTIALITY (Be Like Country)", size(small)) legend(rows(1) size(small)) xsize(12) ysize(11) fysize(60) name(scat1)
logit pref_nuke_A_bin treatment_A##approve_deathpen male age_over60 relig_imp_bin college_grad support_conservative [pweight= weight] if experiment==1 & country!=5 & treatment_A!=5
set scheme s2mono
margins treatment_A#approve_deathpen
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.2)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% prefer nuclear") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title(" " "B. DEATH PENALTY", size(small)) legend(rows(1) size(small)) xsize(12) ysize(11) fysize(51) name(scat2)
graph combine scat1 scat2, imargin(medium) col(1)


**** interactions
* compatriot partiality (own vs libya)
logit pref_nuke_A_bin treatment_A##own_vs_Libya_feeling_50plus male age_over60 relig_imp_bin college_grad support_conservative [pweight= weight_new_IL_therm_any] if country!=5 & (treatment_A==1 | treatment_A==2)
margins treatment_A#own_vs_Libya_feeling_50plus
marginsplot, plotregion(margin(8 8 2 2)) plotopts(mlabel(_margin) mlabangle(45) yvarformat(%5.4f) mlabsize(vsmall)) title("preference for nuclear strike in conditions 1 and 2, by compatriot preference", size(small)) xtitle("treatment group") ytitle("% prefer nuclear") xlabel(.8 " " 1 "effectiveness 90/90" 2 "effectiveness 90/45" 2.20 " ", notick labsize(small)) legend(col(1))
contrast treatment_A#rb1.own_vs_Libya_feeling_50plus, effects

* compatriot partiality (own vs libya) with CONTINUOUS feeling thermometer
logit pref_nuke_A_bin treatment_A##c.own_vs_Libya_feeling male age_over60 relig_imp_bin college_grad support_conservative [pweight= weight_new_IL_therm_any] if (treatment_A==1 | treatment_A==2) & country!=5
margins treatment_A, at(own_vs_Libya_feeling=(-100(10)100)) vsquish
marginsplot, recast(line) recastci(rarea) yline(0) plotopts(mlabel(_margin) mlabangle(45) ylabel(#10,grid labsize(small)) yscale(range(0 .4)) yvarformat(%5.2f) mlabsize(vsmall)) ytitle("effect on probability of preferring the nuclear strike") xtitle("Own minus Libya feeling thermometer") xlabel(, labsize(vsmall)) title("Effect of Own vs. Libya feeling on preference for nuclear weapons" "from A1 (90/90) to A2 (90/45)", size(3))xsize(12) ysize(8)

* death penalty
logit pref_nuke_A_bin treatment_A##approve_deathpen male age_over60 relig_imp_bin college_grad support_conservative [pweight= weight] if country!=5 & (treatment_A==1 | treatment_A==2)
margins treatment_A#approve_deathpen
marginsplot, plotregion(margin(8 8 2 2)) plotopts(mlabel(_margin) mlabangle(45) yvarformat(%5.4f) mlabsize(vsmall)) title("preference for nuclear strike in conditions 1 and 2, by death penalty approval", size(small)) xtitle("treatment group") ytitle("% prefer nuclear") xlabel(.8 " " 1 "effectiveness 90/90" 2 "effectiveness 90/45" 2.20 " ", notick labsize(small)) legend(col(1))
contrast treatment_A#rb1.approve_deathpen, effects

** distribution of compatriot partiality
* in conditions 1-4
mean own_vs_Libya_feeling_50plus, over( country_ascending), [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 & treatment_A!=5
* in control only
mean own_vs_Libya_feeling_50plus, over( country_ascending), [pweight= weight_new_IL_therm_any] if treatment_C==1 & country!=5

* distribution of death penalty
* in conditions 1-4
mean approve_deathpen, over( country_ascending), [pweight= weight] if experiment==1 & country!=5 & treatment_A!=5
* in control only
mean approve_deathpen, over( country_ascending), [pweight= weight] if treatment_C==1 & country!=5


*************************** Figure 5 Marginals for nuclear preference and ethics
drop _est_prefer 
drop _est_ethics
logit pref_nuke_A_bin france israel UK USA male age_over60 relig_imp_bin college_grad support_conservative approve_deathpen own_vs_Libya_feeling_50plus treatmentA_2 treatmentA_3 treatmentA_4 treatmentA_1 [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 & treatmentA_5==0
margins, dydx(*) post
est store prefer
logit nuke_ethical_A_bin france israel UK USA male age_over60 relig_imp_bin college_grad support_conservative approve_deathpen own_vs_Libya_feeling_50plus treatmentA_2 treatmentA_3 treatmentA_4 treatmentA_1 [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 & treatmentA_5==0
margins, dydx(*) post
est store ethics
coefplot (prefer, msymbol (S)) (ethics, msymbol (D)), omitted xline(0) msymbol(.) ciopts(recast(rcap)) yvarformat(%0.1f) yscale(range(0 1)) ytitle("variable", height(5)) order (UK USA france israel male age_over60 relig_imp_bin college_grad support_conservative approve_deathpen own_vs_Libya_feeling_50plus treatmentA_1 treatmentA_2 treatmentA_3 treatmentA_4) ylabel(1 "UK" 2"USA" 3 "France" 4 "Israel" 5 "Male" 6 "Age over 60" 7 "Religiosity" 8 "College Grad" 9 "Conservative" 10 "Approve Death Penalty" 11 "High Compatriot Partiality" 12 "Condition 1" 13 "Condition 2" 14 "Condition 3" 15 "Condition 4")  xtitle("marginal effect", height (8)) ysize(6) xsize(4)


*** ADDITIONAL EQUATIONS

** ALTERNATE FIGURE 2 (including Christian condition)
graph drop _all
logit pref_nuke_A_bin treatment_A##country_ascending [pweight= weight] if experiment==1 
set scheme s2mono
margins, over (treatment_A country_ascending)
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.8)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% prefer nuclear") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title("A. If you had to choose between one of the two military options described in the news article," "would you prefer the nuclear strike or the conventional strike?", size(small)) legend(off) xsize(17) ysize(11) fysize(51) name(scat1)
logit nuke_ethical_A_bin treatment_A##country_ascending [pweight= weight] if experiment==1
set scheme s2mono
margins, over (treatment_A country_ascending)
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.8)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% nuclear ethical") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title("B. Regardless of which option you preferred, how ethical or unethical" "do you think it would be if [country] decided to conduct the nuclear strike in this situation?", size(small)) legend(off) xsize(17) ysize(11) fysize(51) name(scat2)
logit conv_ethical_A_bin treatment_A##country_ascending [pweight= weight] if experiment==1
set scheme s2mono
margins, over (treatment_A country_ascending)
mplotoffset, offset(0.16) recast(scatter) plotregion(margin(5 5 0 0)) plotopts(mlabel(_margin) mlabangle(0) yvarformat(%5.2g) mlabsize(1.8) msize(1.8)) ylabel(#10,grid labsize(small)) yscale(range(0 1)titlegap(*12)) ytitle("% conv. ethical") xtitle(" ") xlabel(, angle(0) alternate labsize(small)) title("C. Regardless of which option you preferred, how ethical or unethical" "do you think it would be if [country] decided to conduct the conventional strike in this situation?", size(small)) legend(rows(1)) xsize(17) ysize(11) fysize(65) name(scat3)
graph combine scat1 scat2 scat3, imargin(vsmall) col(1)

** ALTERNATE FIGURE 5 (including "terrorism likely" control)
drop _est_prefer 
drop _est_ethics
logit pref_nuke_A_bin france israel UK USA male age_over60 relig_imp_bin college_grad support_conservative approve_deathpen own_vs_Libya_feeling_50plus treatmentA_2 treatmentA_3 treatmentA_4 treatmentA_1 terrorism_likely_A_bin [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 & treatmentA_5==0
margins, dydx(*) post
est store prefer
logit nuke_ethical_A_bin france israel UK USA male age_over60 relig_imp_bin college_grad support_conservative approve_deathpen own_vs_Libya_feeling_50plus treatmentA_2 treatmentA_3 treatmentA_4 treatmentA_1 terrorism_likely_A_bin [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 & treatmentA_5==0
margins, dydx(*) post
est store ethics
coefplot (prefer, msymbol (S)) (ethics, msymbol (D)), omitted xline(0) msymbol(.) ciopts(recast(rcap)) yvarformat(%0.1f) yscale(range(0 1)) ytitle("variable", height(5)) order (UK USA france israel male age_over60 relig_imp_bin college_grad support_conservative approve_deathpen own_vs_Libya_feeling_50plus treatmentA_1 treatmentA_2 treatmentA_3 treatmentA_4) ylabel(1 "UK" 2"USA" 3 "France" 4 "Israel" 5 "Male" 6 "Age over 60" 7 "Religiosity" 8 "College Grad" 9 "Conservative" 10 "Approve Death Penalty" 11 "High Compatriot Partiality" 12 "Condition 1" 13 "Condition 2" 14 "Condition 3" 15 "Condition 4" 16 "terrorism likely")  xtitle("marginal effect", height (8)) ysize(6) xsize(4)

**manipulation checks
tab manip1_effective_A treatment_A, col, if treatmentA_5==0
tab manip2_A treatment_A, col, if treatmentA_5==0

**Mass Terrorism Likely
mean terrorism_likely_A_bin, over(country), [pweight= weight] if experiment==1 & country!=5, coefleg
test _b[c.terrorism_likely_A_bin@3.country] = _b[c.terrorism_likely_A_bin@4.country]

logit pref_nuke_A_bin terrorism_likely_A_bin USA france israel UK male age_over60 relig_imp_bin college_grad support_conservative approve_deathpen own_vs_Libya_feeling_50plus treatmentA_5 treatmentA_2 treatmentA_3 treatmentA_4 treatmentA_1 [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 
logit nuke_ethical_A_bin terrorism_likely_A_bin USA france israel UK male age_over60 relig_imp_bin college_grad support_conservative approve_deathpen own_vs_Libya_feeling_50plus treatmentA_5 treatmentA_2 treatmentA_3 treatmentA_4 treatmentA_1 [pweight= weight_new_IL_therm_any] if experiment==1 & country!=5 

*** correlates of feeling thermometer missingness in Israel sample
logit country_feeling_missing pref_nuke_A_bin if country==4
logit country_feeling_missing nuke_ethical_A_bin if country==4
logit libya_feeling_missing pref_nuke_A_bin if country==4
logit libya_feeling_missing nuke_ethical_A_bin if country==4


* terrorism likely by country
mean terrorism_likely_A_bin , over( country_ascending), [pweight= weight] if experiment==1 & country!=5 & treatment_A!=5

* country "must protect its most vital interests even if that sometimes requires doing things that are morally wrong" (conditions 1-4)
mean interests_over_morals_A, over(country_ascending), [pweight= weight] if experiment==1 & treatmentA_5==0

* country "must protect its most vital interests even if that sometimes requires doing things that are morally wrong" (control conditrion only)
mean interests_over_morals_C, over(country_ascending), [pweight= weight] if experiment==3














